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ABSTRACT 


Global soil moisture and land-surface evapotranspiration fields are 
computed using an analysis scheme based on the Simple Biosphere (SiB) soil- 
vegetation-atmosphere interaction model. The scheme is driven with observed 
precipitation, and potential evapotranspiration, where the potential 
evapotranspiration is computed following the surface air temperature- 
potential evapotranspiration regression of Thornthwaite (1948). The observed 
surface air temperature is corrected to reflect potential (zero soil moisture 
stress) conditions by letting the ratio of actual transpiration to potential 
transpiration be a function of normalized difference vegetation index (NDVI), 
based on advanced very high resolution radiometer (AVHRR) data. The land- 
surface hydrology model features the integration of soil moisture equations 
and simplified evaporation routines patterned after SiB. Soil moisture, 
evapotranspiration, and runoff data are generated on a daily basis for a 10-year 
period, January 1979 through December 1988, using observed average monthly 
precipitation scaled by the frequency of observed 1979 daily precipitation events, 
and gridded at a 4° by 5° resolution. The data set provides model-compatible 
soil moisture initial conditions for general circulation models (GCMs) using 
the SiB land-surface parameterization. In addition, the methodology described 
can be readily adapted to compute similarly compatible soil moisture 
initialization fields for other land-surface schemes. 
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1. Introduction 


Several numerical studies have shown that large-scale atmospheric 
circulation and precipitation are sensitive to land-surface evapotranspiration 
(e.g., Shukla and Mintz, 1982; Rowntree, 1983; Yeh et al., 1984; Sud and Smith, 
1985; and others reviewed by Mintz, 1984). In addition, it is now generally 
recognized that land-surface evapotranspiration is one of the more important 
physical processes of the global weather and climate system - and that soil 
moisture, which strongly influences evapotranspiration, is an important 
physical state variable of the system. Soil moisture and land-surface 
evapotranspiration can be determined by employing a moisture budget model 
and knowledge of the potential evapotranspiration (the evapotranspiration in 
the absence of water stress) obtained by applying the conservation of energy 
equation at the Earth’s surface. 

Penman (1948), Budyko (1956), and Priestley and Taylor (1972) showed 
that the potential evapotranspiration is most strongly dependent upon the net 
all- wavelength downward radiation flux, R„, at the surface of the Earth. This 
radiation flux is a function of cloud structure, surface albedo, and the vertical 
distribution of temperature and water vapor in the atmosphere. The vertical 
distributions of atmospheric temperature and moisture can be estimated from 
satellite retrievals and radiosonde measurements. Unfortunately, however, 
the rapidly varying cloud structure and land-surface temperature introduce 
complexities that make an accurate determination of the surface radiation 
balance difficult. In an effort to circumvent some of these radiation-related 
difficulties, Thornthwaite (1948) introduced a method to determine potential 
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evapotranspiration using surface air temperature. This method invokes an 
empirical regression between the potential evapotranspiration, the observed 
daily mean surface air temperature over well-watered vegetation, and the 
length of daylight; employing the hypothesis that the air temperature under 
potential conditions and the daylight length serve as a proxy measurement for 
R n (Thomthwaite and Hare, 1965). Willmott et al. (1985) showed that the 
Thornthwaite procedure produces monthly potential evapotranspiration 
estimates with root mean square errors that are comparable to those produced 
by several other commonly applied radiation-based and combination methods 
(Jensen, 1973; Willmott, 1984). In a recent paper, Mintz and Walker (1993) 
compared the potential evapotranspiration computed using Thomthwaite’s 
methodology with the energy balance schemes of Penman, Budyko, and 
Priestly-Taylor at locations where net radiation flux measurements were 
available. They concluded that potential evapotranspiration computed using 
the Thornthwaite equation closely agrees with that computed by the net 
radiation schemes (to within ± 0.5 mm day 1 ). They also confirmed that 
Thornthwaite’s method provides a reasonable representation of 
evapotranspiration under potential conditions. In view of these findings, we 
use the Thomthwaite (1948) regression analysis to determine global fields of 
potential evapotranspiration. 

In the Simple Biosphere Model (SiB) of Sellers et al. (1986) and several 
other land-surface hydrology parameterizations (e.g., Abramopoulos et al., 
1988; Wood et al., 1992), computed evapotranspiration, and subsequently all 
other components of the surface energy balance, are strongly dependent upon 
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soil moisture. Consequently, the initialization of soil moisture plays a key role 
in the evolution of simulated fields during the early phases of a general 
circulation model (GCM) integration (e.g., Sato et al., 1989a). Serafini and Sud 
(1987) show soil moisture time scales ranging from weeks to several months (5 - 
90 days), depending upon the season and geographic region. Global soil 
moisture fields used for initialization of GCMs are generally produced by 
simple water budget calculations (e.g., Willmott, 1977). In these schemes, the 
actual evapotranspiration is typically defined to be proportional to the potential 
evapotranspiration, where the proportionality factor, or P function, is 
dependent upon soil moisture. Fennessy and Sud (1983), for example, defined 
P to be the soil moisture function based on Nappo (1975). In their study and 
several others (e.g., Willmott et al., 1985; Schemm et al., 1992; Mintz and 
Walker, 1993), the potential evapotranspiration is computed based on 
Thomthwaite (1948). 

The complexities of biophysical control over moisture exchanges 
between the soil and atmosphere in the Simple Biosphere Model suggests that 
an internally consistent soil moisture initialization procedure will include 
much of the same control over these processes as the actual model. Sato et al. 
(1989b) described a method to develop SiB-compatible soil moisture that 
transforms soil moisture fields, such as those produced by Willmott et al. 
(1985) and Gordon (1986), into SiB-compatible fields by preserving some of the 
SiB vegetation control over evapotranspiration; however, they did not take into 
consideration the influence of the moisture transport within the modeled soil 
column. Thus, a global soil moisture analysis scheme that is fully consistent 
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with the soil hydrology of SiB is still needed to produce SiB- compatible soil 
moisture initialization fields. Such a scheme would, at a minimum, feature 
the integration of SiB soil moisture and evapotranspiration equations with 
observed precipitation and estimated potential evapotranspiration. 

Measured surface air temperature is frequently obtained under ‘non- 
potential’ surface conditions. To compute potential evapotranspiration using 
Thomthwaite (1948), the observed surface air temperature requires correction 
to reflect potential (zero soil moisture stress) conditions; if such a correction is 
not made, and observed surface air temperatures are used directly in the 
Thomthwaite equation, it can lead to a misrepresentation of the potential 
evapotranspiration, in part because the reduction in evaporation leads to 
increased sensible heat and, consequently, higher air temperatures than those 
found under potential conditions. In our formulation, this temperature 
correction is accomplished by letting the ratio of actual transpiration to 
potential transpiration be a function of normalized difference vegetation index 
(NDVI), determined from satellite measurements. This eliminates the need to 
consider root zone soil moisture and the soil moisture storage capacity when 
performing the temperature correction, as was done, for example, by Mintz 
and Walker (1993). Daily, 4° by 5° gridded, global evapotranspiration, runoff, 
and soil moisture for each of the three SiB soil layers, are generated for the 10- 
year period January 1979 through December 1988, using observed fields of 
surface air temperature, precipitation, and NDVI. 
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2. Land-surface hydrology model 

A simple local water budget for the land surface can be written as 

^=P + I- E- Q (1) 

at 

where S is the storage of surface and groundwater, P is the precipitation rate, 
E is the evapotranspiration rate, I is the inflow rate from adjacent regions, Q is 
the runoff rate or discharge from surface and groundwater, and t is time. For 
modeling purposes, the storage term in (1) can be cast in the form 

S=W0 8 D (2) 

where W is the fractional soil wetness, 0 < W < 1, 0, is the porosity (or the 
volumetric soil moisture at saturation), and D is the hydrologically active soil 
column depth. 


a. Soil moisture 

The foundation of the soil moisture computations performed in this 
study is the discretization of (1) and (2) into the three-layer soil moisture 
equations used in SiB. The land area contained within each GCM grid box is 
considered horizontally uniform, but the soil— vegetation types vary globally. 
In this application, the SiB governing soil moisture equations have been cast in 
the form 


aw, = 1 
at QgD, 
aw 2 i 
at e 8 D 2 
dW 3 _ 1 

at e 8 D 3 


[P> - Q« - E,] 

[Qi,2 — Q2,3 — ^ 2 ] 
[Q2.3 ~ Q3 ~ Qb] 


(3a) 

(3b) 

(3c) 
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where the subscripts 1, 2, and 3 indicate the three soil layers; Q i(i+l is the 
moisture transfer rate between layers i and i+1, for i = 1, 2; Q 3 is the 
gravitational drainage rate from layer three; E x is the soil evaporation rate; E 2 
is the rate of root zone moisture loss due to transpiration; and Q b is baseflow 
drainage rate produced by a linear reservoir model. These equations are 
similar to those presented by Sellers et al. (1986) except that they contain fewer 
evaporation terms in layers one and two and have an additional moisture 
drainage term, Q b , in layer three. 

The general flux-moisture gradient structure of this three-layer soil 
hydrology model is similar to several other state-of-the-art land- surface 
hydrology parameterizations used in GCMs, including the Biosphere- 
Atmosphere Transfer Scheme (BATS) (Dickinson et al., 1986), the Goddard 
Institute for Space Studies model (Abramopoulos et al., 1988), and the 
Canadian Land Surface Scheme (CLASS) (Verseghy, 1991). The moisture 
exchange between soil layers is determined by the steady-state, unsaturated, 
one-dimensional solution of Darcy’s law (Freeze and Cherry, 1979), 

« ' -«(g * ■) 

where Q is the downward transport of water, K is the hydraulic conductivity of 
the soil, v is the soil moisture potential, and z is the vertical coordinate. To 
solve for the soil moisture transport, (4) is discretized, and the soil moisture 
potential, \|/i> of layer i, is given by 

v, = v. w,- B (5) 

where \y B is the soil moisture potential at saturation, and the coefficient B is a 
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function of soil type (Clapp and Homberger, 1978). The hydraulic conductivity, 
Kj, of layer i, is given by (Campbell, 1974) as 

K, = K 8 W/ 2B+3) (6) 

where the saturated hydraulic conductivity, K*, varies with soil characteristics 
as listed by Clapp and Homberger (1978). Gravitational drainage transports 
moisture out of the lowest layer, such that 

Q 3 = -K 3 sin(p) (7) 

where p is the mean slope angle. 

In addition to the gravitational drainage from the third soil layer, 
baseflow runoff from layer three is produced by a linear reservoir model. In 
this formulation, subgrid-scale spatial variability of soil moisture has been 
included in the runoff formulation by assuming that the grid box contains a 
saturated areal fraction that is equal to the soil wetness of the third soil layer 
(Liston et al., 1993). This saturated region then drains as a linear reservoir, 
thus producing a term, Q b , in the governing equation for the third soil 
moisture store, 

Q b = k b W 3 (8) 

where k b is the discharge rate from the third layer at saturation. In SiB, there 
are twelve soil-vegetation categories corresponding to the indices 1 through 12 
(Dorman and Sellers, 1989). Based on off-line simulations over the Mississippi 
River Basin, for ‘grassland’ (soil-vegetation type 7), k b 7 was shown to be 0.6 mm 

day' 1 (Liston et al., 1993), where the number superscript indicates vegetation 
type index. To allow extrapolation to other soil-vegetation types, k b is scaled by 
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the depth of layer three, D 3 \ for the j th soil-vegetation type, according to 



where 1 < j < 12. 

In this implementation, all precipitation reaching the ground enters the 
soil matrix. In the original SiB model, the precipitation was given only one 
time step to infiltrate the soil. Any moisture that did not infiltrate during that 
time step, called precipitation excess, was identified as surface runoff. This 
was found to produce unrealistically large runoff events during periods of low 
hydraulic conductivity, such as late summer, when the precipitation could not 
infiltrate the soil. Using a runoff routing model designed to generate river 
discharge hydrographs from gridded runoff data produced by a land-surface 
hydrology model or GCM (Liston et al., 1993), these Hortonian-like runoff 
events, occurring at the 4° by 5° GCM scale, were so significant that they 
produced an anomalous autumn peak in the Mississippi River Basin 
discharge hydrograph. 

In summary, the significant changes made to the original SiB soil 
hydrology formulation are first, addition of linear reservoir drainage, through 
the inclusion of the drainage term Qb in (3c), and its description given by (8) 
and (9); and second, the practice of forcing all precipitation to enter the upper 
soil matrix in one time step. These modifications to the SiB soil hydrology are 
currently being tested in on-line simulations using the Goddard Laboratory for 
Atmospheres (GLA) SiB-GCM. 
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b. Evapotranspiration 

Soil moisture dependent evapotranspiration is computed by scaling the 
potential evapotranspiration, E p , by a factor ft, 

E s = ft E p (10) 

where the subscript i = 1, 2 identifies the soil layer of interest, and E is the 
actual evaporation rate. For the case of vegetated land, ft is defined to be the 
soil water stress factor for the second soil layer, flty 2 ), given by Xue et al. (1991), 

f(V 2 ) = 1 - exp(-c 2 (cj - ln(-\jr 2 ))) 0<f(\|/ 2 )<l (11) 

where c x and c 2 are empirical constants representing the progression of 
closing stomata. Values of Ci and c 2 for the different vegetation types are given 
in Table 1, where the swapping of Ci and in the table presented by Xue et al. 
(1991) has been corrected. For the case of soil-vegetation type 11 (bare soil), ft = 
Wj, the soil wetness of the uppermost soil layer. 

c. Atmospheric forcing 

In its reduced form described above, the SiB land-surface hydrology 
parameterization was forced using daily values of potential evaporation and 
surface station precipitation climatology. Potential evapotranspiration, E p , 
was computed following Thorn thwaite (1948), which requires a surface air 
temperature representative of potential conditions. To account for measured 
surface air temperatures, that may not have been collected under conditions of 
zero soil moisture stress, we let the ratio of actual transpiration to potential 
transpiration be a function of normalized difference vegetation index (NDVI), 
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available from satellite measurements. This has eliminated the need to 
consider root zone soil moisture and the soil water storage capacity when 
computing the temperature correction, and has also prevented errors in 
measured precipitation from affecting the computations (Mintz and Walker, 


1993). 

In Thornthwaite’s empirical formulation, the daily potential 
evapotranspiration is a function of the number of daylight hours, h, and the 


mean daily air temperature, T,\ obtained under conditions of no water stress, 


E 


p 


0 , 


T a * < 0°C 



(-13.86 + 1075 T a * - 0.0144 T a * 2 )(% 2 ) 


0<T a * < 26.5°C 
26.5°C<T a * (12) 


with 



T m > 0°C (13) 


and 

a = (6.75 xlO - 7 I 3 ) - (7.71X10" 6 I 2 ) + (179X10" 2 i) + 0.492 (14) 
where T m is the monthly mean surface air temperature; and I is the annual 
heat index, or sum of the 12 monthly heat indices, i. The daylight duration, in 


hours, is calculated from 

h = (24/180) arc cos(- tan<J> tanS) 
where <J> is latitude, and 5 is the Sun’s declination, given by 


8 — (J)*p cos 


/ 

! 

1 

T3 

2k | 

d II 

V 

V. y JJ 


(15) 


(16) 
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where is the latitude of the Tropic of Cancer (23.45°), d is the Gregorian Day, 
d r is the day of the summer solstice (173), and d y is the average number of days 
in a year (365.25). Poleward of 50° latitude, Thomthwaite defined h to be the 
same as h at 50°. 

The measured surface air temperature T a , is related to the air 
temperature under potential conditions according to 

T a * = (0.83 + 0.17 J3 e )T a (17) 

where (3 e is a soil moisture dependent evapotranspiration coefficient. This 
relationship was derived from a comparison of the measured daily mean 
surface air temperatures in moist humid regions, and in dry desert regions at 
the same latitude and the same season of the year (Mintz and Walker, 1993). 

In describing the evapotranspiration coefficient, p e , Mintz and Walker 
considered the evaporative balance, 

E t + E s = P e (E p - E,) (18) 

where E T , E s , and E r are the transpiration, soil evaporation, and interception 
losses, respectively. They defined p e to be a function of soil wetness, 

p e = 1 - exp (- 6.8 (W/W*)) (19) 

where W is the root zone soil moisture, and W* is the root zone water storage 
capacity (Davies and Allen, 1973; Nappo, 1975). 

In our application we follow Mintz (unpublished manuscript) and 
combine the interception loss (precipitation accumulated on the vegetative 
canopy and evaporated) and the soil evaporation into a single component, Ei >s , 
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assuming that they behave in a roughly complementary fashion. When 
precipitation is greater than the potential evapotranspiration, increasing 
(decreasing) leaf area index will increase (decrease) the amount of 
precipitation intercepted by, and evaporated from, the leaves, and there will be 
a corresponding reduction (increase) in the amount of insolation reaching the 
soil for soil evaporation. Thus, we let 

fl, = E, / (E, - E, s ) (2° 

where Ei s is considered to be the lesser of the daily precipitation, P , and daily 
potential evapotranspiration. 

To derive an equation for the evapotranspiration coefficient, assume that 
E t / (E p - Eis) is proportional to the area average stomatal conductance, l/r g , or 
the reciprocal of the stomatal resistance, r g , 


E r /(E, - E,, s ) - 1/r. (21) 

Adopting the findings of Sellers (1985), we let the area average stomatal 
conductance be proportional to the simple ratio, In<NiR> I lH(vis)> 

^l(NIR) 


1/r. 


In 


( 22 ) 


l(VIS) 


where In (N iR) and Inorcs) are the nadir or near-nadir surface reflectance in the 
near infrared (0.7-3.0 pm) and visible wavelenth (0.4-0.7 pm) intervals, 
respectively. This simple ratio is related to the normalized difference 
vegetation index, NDVI (Tucker et al., 1985), according to 


Iri(NiR) _ (1 + NDVI) _ 
lTl ( vis) " (1 - NDVI) ’ 


NDVI 


_ (^(NIR) IH(VIS)) 
(l^NIR) “ ^(VIS)) 


(23) 
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Combining (20) - (23) leads to a formulation for the evapotranspiration 
coefficient (3 e , which is a function of NDVI, 


Pe = Y + * 


(1 + NDVI) 


(24) 


(1 - NDVI) 

where the coefficients y and X can be obtained by letting p e = 0 and 1 when 
NDVI equals the lower and upper limits of photosynthetically active 
transpiring vegetation, respectively. For example, values of y = -0.157 and X = 
0.129 result if a lower NDVI limit of 0.10 and an upper limit of 0.80 are 
assumed. 

It would be desirable to use daily precipitation observations in the model 
integrations. Unfortunately, the unavailability of global daily precipitation 
data for the 10-year period of interest precludes such an implementation. An 
alternative is to use the available daily global precipitation record from the 1979 
First GARP Global Experiment (FGGE) year as a stencil for the other years of 
the period. Monthly station precipitation climatology (Spangler and Jenne, 
1990) for the 10-year period, 1979 through 1988, were reduced to a 4° latitude by 
5° longitude grid. A surrogate temporal record of daily precipitation was 
produced by scaling the gridded monthly precipitation climatology, by the 
magnitudes of the 1979 FGGE daily precipitation data, using the formula 


{ x> c \ 


P - P 

r d - r d 


P F 

V r m J 


(25) 


where P d c is the daily precipitation, P d F is the observed daily precipitation in 
the FGGE year, P m F is the observed monthly precipitation in the FGGE year, 
and P m c is the observed monthly precipitation climatology. Thus, the daily 
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precipitation forcing preserves measured monthly magnitudes and occurs at 
the frequency of the observed 1979 precipitation events. Admittedly, in the 
natural system, the frequency of precipitation events will vary from one year to 
the next. In this application, we have assumed that when the daily computed 
soil moisture, evapotranspiration, and runoff are compiled to produce monthly 
means, the resulting values will be representative of the period in question. 
The availability of daily precipitation data over the United States is currently 
being used to test this assumption. In the integrations, all precipitation data 
were treated as rainfall. Station monthly mean surface air temperature data 
were gridded and linearly interpolated to provide mean daily air temperatures. 

3. Presentation of data 

By aggregating the daily data produced by the model integrations, global 
monthly mean precipitation, evapotranspiration, runoff, and soil moisture for 
each of SiB’s three soil layers have been produced for each month of the 10-year 
period, 1979 through 1988. In this technical memorandum, global distribution 
maps illustrating 10 -year (1979 - 1988) monthly averages of the monthly mean 
precipitation, evapotranspiration, root zone soil moisture, and runoff fields 
have been plotted. In addition, the standard deviations of the monthly mean 
precipitation, evapotranspiration, root zone soil moisture, and runoff for the 10- 
year period are also included. These are found in the section titled: 
“HYDROLOGIC FIELDS: 10-Year Average of Monthly Means, and Standard 
Deviations” under the subheadings, PRECIPITATION (January - December), 
EVAPOTRANSPIRATION (January - December), ROOT ZONE SOIL 
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WETNESS (January - December), and RUNOFF (January - December). 

Referring to the precipitation figures, in January, maximum 
precipitation values of over 8 mm day' 1 occur in the tropical regions south of 
the equator in South America, Africa, the Pacific Islands, and Australia. In 
addition, a January precipitation maximum is found in the mountainous 
regions surrounding the Gulf of Alaska. In January, regions of precipitation 
less than 0.5 mm day 1 include Arctic North America, the central United 
States, central Africa, and central Asia. By July, the tropical precipitation 
peak has shifted north of the equator to include South and Central America, 
Africa, and India. The precipitation minimum in central Africa has shifted 
north. A minimum of less than 0.5 mm day' 1 in July also exists in southern 
Africa, central South America, and central Australia. In general, the 
precipitation standard deviation patterns correspond to the precipitation 
patterns, with the highest standard deviations occurring in the regions of 
highest precipitation. 

Evapotranspiration is minimal, less than 0.5 mm day' 1 , in the Northern 
Hemisphere in January, while south of the equator, evapotranspiration is 
typically greater than 2 mm day' 1 . In May, regions north of 30° N latitude 
show a significant increase in evapotranspiration. By July, Northern 
Hemisphere evapotranspiration has increased sharply, typically being 
between 2 and 4 mm day ’ 1 , with the exception of the regions of northern Africa 
and Arabia, where values less than 0.5 mm day' 1 are typical. In general, the 
evapotranspiration standard deviation patterns correspond to the 
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evapotranspiration patterns, with, the highest standard deviations occurring 
in the regions of highest evapotranspiration. Typically the evapotranspiration 
standard deviation displays more structure than the relatively smooth 
evapotranspiration patterns. In the Southern Hemisphere winter, South 
America has minimal evapotranspiration standard deviation. 

Referring to the soil wetness figures, the lowest January root zone soil 
wetness values are found in the southwest United States, Mexico, southern 
South America and Africa, the Sahara, Arabia, and central Asia. These 
minimum soil moisture fractions are typically between 0.15 and 0.2. 
Maximum values, between 0.6 and 0.8, are found in the tropics, in the middle 
latitude storm belt regions such as the northwest coastal regions of the United 
States and Canada, and in Europe. In July, subtle changes in the soil 
moisture patterns are evident over the globe, but nothing as dramatic as the 
differences found in the precipitation and evapotranspiration patterns. In 
general, the soil wetness standard deviations vary considerably. Note that in 
the SiB evapotranspiration formulation, although it does depend on soil- 
vegetation type, evapotranspiration occurs at the potential rate for soil wetness 
values around 0.60 and evapotranspiration becomes zero at a soil wetness 
value of near 0.30. As a consequence, a 0.10 change in soil wetness is able to 
have a significant influence on evapotranspiration. Typically the soil wetness 
standard deviations are greatest in the tropics and decrease with increasing 
latitude. 

Peak January runoff occurs in the region surrounding the Gulf of 
Alaska and in the tropics of South America, Africa, the Pacific Islands, and 
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central Australia. In July, the tropical runoff maximum has shifted north, 
and India and regions to the north and east are producing runoff values 
greater than 8 mm day 1 . Since runoff is strongly a function of soil moisture, 
the relatively subtle seasonal changes in soil moisture are reflected in the 
subtle seasonal changes in runoff. The one significant exception to this is the 
dramatic seasonal change in runoff produced by the increased precipitation 
during the Indian monsoon. In general, the runoff standard deviation 
patterns correspond to the runoff patterns, with the highest standard 
deviations occurring in the regions of highest runoff. Compared to the 
minimal runoff standard deviations occurring in the Northern Hemisphere 
during the Northern Hemisphere summer, the tropics contain relatively high 
standard deviations during this period. 

As part of the model integrations, global monthly mean soil moisture for 
each of SiB’s three soil layers were produced for each month of the 10 -year 
period. These soil moisture fields are currently being used as soil moisture 
initial conditions during climate simulations using the version of the GLA 
GCM that employs the SiB land-surface parameterization. 

4. Concluding summary 

An analysis scheme to generate global evapotranspiration and soil 
moisture fields, using key components of the Simple Biosphere Model, has 
been developed. Two modifications have been made to the original SiB soil 
hydrology formulation: a linear reservoir drainage term has been added to the 
third soil layer, and all precipitation is forced to enter the upper soil matrix 
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during each time step. The influence of these modifications is being tested in 
the Goddard Laboratory for Atmospheres SiB-GCM. In addition, the analysis 
scheme incorporates a simplified computation of evapotranspiration based on 
SiB evapotranspiration functions. 

The method requires atmospheric forcing of precipitation and potential 
evapotranspiration. Potential evapotranspiration is computed using the 
method of Thomthwaite, where the observed air temperatures have been 
adjusted to represent potential conditions by using normalized difference 
vegetation index (NDVI) data determined from satellite observations. The 
analysis was run for 10 years, 1979 through 1988. In performing the model 
integrations, it would be desirable to use global observed daily precipitation. To 
circumvent the unavailability of global daily precipitation data for the period 
1979 through 1988, we use the daily global precipitation record from the 1979 
FGGE year as a stencil for the other years of the period. Thus, by using 
observed monthly mean precipitation data, we produce a 10-year, daily 
precipitation record which preserves measured monthly magnitudes, and 
occurs at the frequency of the observed 1979 FGGE year precipitation events. 

The scheme represents an efficient methodology for generating soil 
moisture initialization fields for GCMs; it can also be applied to other land- 
surface hydrology parameterizations being used in GCMs, through a 
substitution of the parameterizations’ soil hydrology and evapotranspiration 
functions. The resulting soil moisture fields are compatible with the 
atmospheric forcing of observed precipitation and surface air temperature 
climatology. In addition, these fields are computed in a manor consistent with 
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the soil hydrology and evapotranspiration representations found in the specific 
GCM land-surface hydrology parameterization of interest. 

Validation of data sets such as the ones described here is an important 
and necessary step in the development and refinement of information used in 
climate simulations using GCMs. At seasonal and continental scales, 
terrestrial precipitation inputs are considered relatively well known, despite 
the scarcity of data in some regions of the world. In addition, river discharge 
measurements are readily available and are being used for model validation 
over large areas (Liston et al., 1993). This river discharge is an important 
integrator of the hydrologic cycle, and is able to provide much needed insight 
into whether or not the models produce realistic terrestrial water balances. 
Evapotranspiration comprises a significant fraction of the moisture transport 
within the hydrologic cycle; unfortunately, direct validation of land- surface 
evapotranspiration is difficult, and it is generally approximated using indirect 
methods that utilize air temperature or radiation observations. Soil moisture 
data from remote sensing measurements are now becoming available for 
comparison with modeled fields. A useful feature of these observations is the 
measurement of subgrid-scale soil moisture variability. These data and other 
satellite-derived products are expected to assist in the formulation of next 
generation runoff and evapotranspiration schemes, and thus provide a 
valuable contribution to future land-surface modeling efforts. The resulting 
improved understanding and description of terrestrial hydrology will add to 
the quality of global climate simulations. 
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NUMERICAL DATA. Input and derived fields described in this 
technical memorandum are available for all months of the 10-year period, 1979 
through 1988, in the form of grid-area numerical values. They are available in 
printed tables, on magnetic tape, and via Internet using standard file transfer 
protocal (ftp). Contact the authors, Climate and Radiation Branch, Code 913, 
Laboratory for Atmospheres, NASA/Goddard Space Flight Center, Greenbelt, 
Maryland 20771. 
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Table 1. Values of the Evapotranspiration Coefficients and Ca for 
the Twelve SiB Vegetation Types* 


Type 

SiB Biome 

Ci 

C2 

l 

Broadleaf-evergreen trees (tropical forest) 

6.25 

1.2 

2 

Broadleaf-deciduous trees 

5.57 

5.35 

3 

Broadleaf and needleleaf trees (mixed forest) 

5.73 

1.92 

4 

Needleleaf-evergreen trees 

5.53 

3.7 

5 

Needleleaf-deciduous trees (larch) 

5.66 

7.8 

6 

Broadleaf trees with grass and shrubs (savanna) 

5.67 

1.8 

7 

Grassland (perennial) 

5.80 

1.73 

8 

Broadleaf shrubs with perennial groundcover 

5.98 

3.0 

9 

Broadleaf shrubs with bare soil (semi-desert) 

6.37 

1.39 

10 

Dwarf trees and shrubs with groundcover (tundra) 

5.37 

0.96 

11 

Bare soil (desert) 

— 

— 

12 

Cultivated land idealized to grow wheat 

4.36 

0.58 


from Xue et al. (1991) 
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HYDROLOGIC FIELDS: 10-Year Average of Monthly Means, 
and Standard Deviations 


Global distribution maps of key hydrologic fields are included in this 
report. Ten-year (1979 - 1988) averages of the monthly mean precipitation, 
evapotranspiration, root zone soil moisture, and runoff fields are plotted. In 
addition, the standard deviations of the monthly mean precipitation, 
evapotranspiration, root zone soil moisture, and runoff for the 10 -year period 
have been included. 


25 




PRECIPITATION 


(January - December) 


Ten-year average (1979 - 1988) of the monthly mean precipitation (mm 
day 1 ), obtained from rain gage observations (Spangler and Jenne, 1990). Also 
plotted are the standard deviations of the monthly mean fields as determined 
from the 10-year data set. 
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Precipitation Standard Deviation (mm/day) 

September 1979-88 
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EVAPOTRANSPIRATION 


(January - December) 


Ten-year average (1979 - 1988) of the monthly mean evapotranspiration 
(mm day' 1 ), computed from soil moisture initialization scheme. Also plotted 
are the standard deviations of the monthly mean fields as determined from the 
10-year data set. 
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EvaDotranspiration Standard Deviation (mm/day) 
K October 1979-88 
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ROOT ZONE SOIL WETNESS 
(January - December) 


Ten-year average (1979 - 1988) of the monthly mean root zone soil 
wetness (0 - 1), computed from soil moisture initialization scheme. Also 
plotted are the standard deviations of the monthly mean fields as determined 
from the 10-year data set. 
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RUNOFF 


(January - December) 


Ten-year average (1979 - 1988) of the monthly mean runoff (mm day' 1 ), 
computed from soil moisture initialization scheme. Also plotted are the 
standard deviations of the monthly mean fields as determined from the 10 -year 
data set. 
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